Neutrino and antineutrino quasielastic interactions with nuclei 
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We investigate the interaction of neutrinos and antineutrinos with nuclei. We explore in particular 
the role played by the multinucleon excitations which can contaminate the quasielastic cross section. 
For neutrinos the multinucleon term produces a sizable increase of the quasielastic cross section. 
Part of the effect arises from tensor correlations. For antineutrinos this influence is smaller owing 
to the axial-vector interference which increases the relative importance of the terms which are not 
affected by these multinucleon excitations. 
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PACS numbers: 25.30.Pt, 13.15.+g, 24.10.Cn 

The progresses of neutrino experimental physics have 
allowed the measurement of several partial cross sections 
in the interaction with nuclei. Most data concern ra- 
tios of cross sections [3-0] but some absolute values are 
now available |8Hl0j. In particular the MiniBooNE col- 
laboration has measured the neutrinos charged current 
(CC) quasielastic (QE) cross section on 12 C for a neutrino 
beam with average energy of 788 MeV [13] ■ In this work 
ejected nucleons are not detected and the quasielastic 
cross section is defined as the one for processes in which 
only a muon is detected in the final state. However it is 
possible that in the neutrino interaction a pion produced 
via the excitation of the A resonance escapes detection, 
for instance because it is reabsorbed in the nucleus, lead- 
ing to multinucleon emission. In this case it simulates 
a quasielastic process. The MiniBooNE analysis of the 
data corrects for this possibility via a Monte-Carlo eval- 
uation of this process. The net effect amounts to a re- 
duction of the observed quasielastic cross section. After 
application of this correction the quasielastic cross sec- 
tion thus defined still displays an anomaly as compared to 
a relativistic Fermi gas prediction. The prediction which 
is sensitive to the cut-off mass of the axial form factor 
fits the data provided a modified axial form factor is in- 
troduced in the calculation, with an increase of the axial 
cut-off mass from the accepted value Ma = 1-03 GeV to 
the value Ma = 1.35_GeV; otherwise the calculated cross 
section is too small (Toj . 

On the theoretical side Martini et al. have drawn 
the attention to the existence of additional sources of 
multinucleon emission which are susceptible to produce 
an apparent increase of the "quasielastic" cross section. 
Their evaluation, although approximate, of this contri- 
bution shows that it is able to account for this apparent 
increase. It stresses in particular the role played by the 
NN tensor correlations in this enhancement. In clarifying 
this point we will be naturally led to explore the corre- 
sponding effect in the case of antineutrinos. We will show 
that antineutrinos can provide a useful test of the origin 
of this anomaly, which is the aim of the present work. 

In our approach, the same as that of Ref. [ll[ , the neu- 
trino cross section on nuclei is expressed in terms of the 
nuclear response functions treated in the random phase 



approximation (RPA) . The only nucleon resonance taken 
into account is the A one. Several responses enter this 
interaction, as exemplified below in a simplified expres- 
sion of the charged current cross section where the lep- 
ton mass is ignored and the A width is taken to zero. It 
reads : 
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where Gf is the weak coupling constant, 9 C the Cabbibo 
angle, k and k' the initial and final lepton momenta, = 
— = (cj, q) the four momentum transferred to the 
nucleus, 9 the scattering angle, M\ (M) the A (nucleon) 
mass. The plus (minus) sign in Eq. (JTJ stands for the 
neutrino (antineutrino) case. The existence of this axial- 
vector interference term is crucial in the present work. 
The various responses are related to the imaginary part 
of the corresponding full polarization propagators : 
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They are related to the inelastic cross section for a given 
coupling. For instance for the isospin operator : 

n j=l 
A 

x(n\J2r(k)e i 1- Xk \0)*S(uj-E n + E ). (3) 
fe=i 

Similar expressions apply to the other transition opera- 
tors. An analogous expression applies for a transition to 
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a state with a A excitation which enters the spin-isospin 
responses. In this case the operators are : 

0* T(i) (j) = (S r q)T± A (T) (j) = (S j xc l ) i T±. (4) 

In the case of a spin operator the index L or T refers 
to the direction of the spin, longitudinal or transverse, 
with respect to the momentum q. In the expression ([1]) 
the upper indices, N or A, refer to the type of particles, 
nucleon or A, excited by the weak current at the two 
ends of the RPA chain (see FigflJ. In the present nu- 
merical evaluations we employ, instead of Eq. ([T]), the 
full expressions given in the Appendix A of our previous 
publication [ll| for which we refer for more details. 

We have treated these responses in the ring approxi- 
mation of the RPA so as to account for collective effects. 
A limitation of our description is that it does not incor- 
porate final state interactions such as for instance the 
possibility for a real pion produced by the neutrino to be 
reabsorbed in the nucleus leading to multinucleon ejec- 
tion. This fraction of the produced pions would thus be 
counted in the "quasielastic" events. These events have 
been subtracted in the MiniBooNE analysis through their 
Monte Carlo evaluation. Hence our theory can be con- 
fronted to their corrected experimental results. 

To lowest order the quasielastic cross section is given 
by the terms in R NN , whether coming from the isovec- 
tor interaction, R^f N , or from the isospin spin-transverse 
one, R^^ T y The isospin spin-longitudinal one is sup- 
pressed, for a vanishing lepton mass, by a cancellation be- 
tween the space and time components of the axial current 
which is the reason why it is not explicitly written 
in Eq. (fTJ) . In the actual calculation its small quasielastic 
contribution is taken into account. It also contribute to 
the multinucleon ejection term. In the RPA chain R NA 
and R AA also contribute to the quasielastic response as 
illustrated in FigfT] The collective effects produce a mild 
suppression of the quasielastic response due to the repul- 
sive nature of the residual interaction. The effect that we 
want to discuss here is of a different nature. It concerns 
the multinucleon (np — nh) ejection in neutrino interac- 
tions other than that due to final state interaction. Here 
several sources of multinucleon emission enter our de- 
scription and all type the responses R NN , R NA and R AA 
contribute. Examples are given in Fig. [5] One of these 
multinucleon sources arises from the modification of the 
A width in the nuclear medium (Fig [2] (d)). This effect 
was introduced and parametrized by Oset and Salcedo 
[l2| in the case of real pion or photon absorption. We 
have used their parametrization of the modified width, 
although the kinematics of neutrino interaction is differ- 
ent since we are in the space-like region, which could be a 
source of uncertainty. For the other terms not reducible 
to a modification of a A width of which some examples 
are shown in Figf5] we have used a parametrization of 
Delorme and Guichon [lijj . They exploited a calculation 
by Shimizu and Faessler [3] of the absorptive part of 
the p-wave pion-nucleus optical potential at threshold, 
which writes (4Tr/2m 7T )V ■ ImCo p 2 V [lj|. It is known 




FIG. 1: Examples of contributions to the quasielastic cross 
section. Lowest order contribution from : (a) R NN ; (b) 
R NA and (c) R AA . The wiggled lines represent the exter- 
nal probe, the full lines correspond to the propagation of a 
nucleon (or a hole), the double lines to the propagation of 
a A and the dashed lines to an effective interaction between 
nucleons and/or As. The dotted lines show which particles 
are placed on-shell. 
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FIG. 2: Examples of contributions to the two nucleon ejection 
term. From : (a) R NN - (b) R NA - (c) and (d) R AA . Diagram 
(d) represents an example of modification of the A width in 
the medium. 



that the absorption mechanism of pions is a two-nucleon 
one which means that in the final state two nucleons are 
ejected. In the many-body language this is a two particle- 
two hole (2p — 2h) excitation. In Ref. the absorption 
is described by three types of terms (see FigJ5]). The first 
one (FigfS] (a) with pion lines replacing the weak current 
ones as in the next diagrams) arises from the nucleon- 
nucleon correlations, essentially from the tensor correla- 
tions. Another one involves a A excitation (Figf2] (c)). 
The third one is an interference between the nucleon cor- 
relation and A terms as in FigJ5] (b). In these graphs 
the coupling of the pion to the nucleon involves the pion 
momentum and is of the spin isospin type, a ■ qf, for the 
nucleon and a similar expression for A excitation. In the 
optical potential the pion momentum is expressed by the 
gradient, the remaining part of the optical potential with 
the parameter ImCo then provides the nuclear 2p — 2h 
bare responses to a probe which couples to the spin and 
isospin of the nucleon. In principle it is a spin longitu- 
dinal coupling but as this term represents a short range 
effect there is no difference in the bare case (i.e. before 
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FIG. 3: Various response contributions to the and z? M - 12 C 
charged current genuine quasielastic cross section. Isovec- 
tor (r) and isospin spin-longitudinal (err (L)) components are 
identical for and z? M . 



the RPA chain) between longitudinal and transverse spin 
couplings. Therefore it applies to the bare magnetic part 
of the vector current and to both the spin-transverse and 
spin- longitudinal parts of the axial current. The pion ab- 
sorption calculation of Ref. [l4| is performed for threshold 
pion, i.e. for a vanishing three-momentum and an en- 
ergy uj = m^, which does not correspond to the neutrino 
situation. Delorme et al. [13j have then introduced in 
each absorption graph the corresponding energy depen- 
dence to obtain the bare 2p — 2h responses to be inserted 
in the RPA chain needed for neutrino interaction. How- 
ever they have completely ignored the momentum depen- 
dence. They have left apart the graph (d) of FigJ3] since 
it corresponds to a modification of the A width which is 
taken into account sep arately through the parametriza- 
tion of Oset et al. [l2[ • In a first approach we have used 
the procedure of Delorme et al. as such in order to evalu- 
ate the bare 2p—2h components to be inserted in the RPA 
chain for the evaluation of the neutrino cross section. In 
our resulting 2p — 2h cross section the modification of 
the A width is not a dominant effect; it adds a small 
2p — 2h component (plus a 3p — 3h one). It turned out 
that our overall multinucleon contribution, when added 
to the genuine quasielastic cross section, was able to ac- 
count for the anomaly without modification of the axial 
cut-off mass as already reported in Ref. [ll| . 

Since our parametrization of the 2p — 2h piece from the 
extrapolation of pion absorption is questionable, as it ig- 
nores in particular any momentum dependence, we have 
also investigated the effect with our second parametriza- 
tion [ll| of the 2p—2h contribution beyond the one which 
is reducible to a modification of the Delta width. In the 
second approach we have used a microscopic calculation 
of Alberico et al. [l6[ specifically aimed at the evaluation 
of the 2p — 2h contribution to the isospin spin-transverse 
response, measured in inclusive (e, e') scattering. Their 
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FIG. 4: Differential CC and - 12 C cross sections versus 
the energy transfer. 



basic graphs are similar to those of Shimizu and Faessler. 
In principle this way of evaluation is definitely more sat- 
isfactory since the kinematical variables are correctly in- 
corporated but the results of Ref.[l6[ are available only 
for a limited set of energy and momenta. We have ex- 
tended this range to cover the neutrino one through an 
approximate extrapolation. We refer to for the de- 
tails. The result was that the distribution in the en- 
ergy transfer w of the differential neutrino cross section 
is largely modified with a more realistic distribution but 
once integrated in energy the difference is small. The cor- 
responding neutrino and antineutrino total "quasielastic" 
cross sections are practically unchanged. 

Now from the way in which this 2p — 2h contribution 
is built, the corresponding coupling of the weak current 
to the nucleon or A is a spin isospin one. It is then clear 
that this 2p — 2h term only affects the magnetic and axial 
responses which enter the neutrino cross section. In the 
expression (fTJ) these are the terms in G\ , G 2 M and the in- 
terference term in GaGm- The isovector response (term 
in R T ) instead is not affected. This difference is the basis 
for the test that we propose to help elucidate the origin of 
the anomaly. The principle is simple. For antineutrinos 
the interference term in GmGa produces a suppression of 
the spin isospin response contribution, while it enhances 
it for neutrinos. In the neutrino case the contribution 
from the spin isospin terms largely dominates the isovec- 
tor one. If the suppression of the spin isospin part for 
antineutrinos by the interference term modifies the bal- 
ance in such a way that the role of the isovector response 
in the cross section becomes appreciable, the relative role 
of the 2p — 2h part will be smaller for antineutrinos. We 
will show below that this is indeed the case. Figure[3]dis- 
plays the various components to the genuine quasielastic 
(single nucleon ejection) cross section on 12 C as a func- 
tion of the energy both for muon neutrinos and antineu- 
trinos. The isovector part is identical in the two cases, 
as well as the isospin spin-longitudinal one. The isospin 
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FIG. 5: (upper panel) and (lower panel) - 12 C CC 
"quasielastic" cross sections per neutron and per proton with 
and without the multinucleon component as a function of neu- 
trino energy. The experimental points are taken from [Tol ]. 



FIG. 6: Ratio of multinucleon component of "quasielastic" 
cross section on 12 C to the single nucleon one for and £> M 
as a function of neutrino energy. 



spin-transverse one instead changes considerably between 
neutrinos and antineutrinos due to the interference term 
with an appreciable reduction in the antineutrino case. 
As a consequence the isovector relative contribution to 
the total cross section becomes quite significant for an- 
tineutrinos which opens the possibility of an experimen- 
tal test. Such details are visible in Fig. 2] where the 
differential cross section 4^ at MeV is shown in 

the bare and RPA case, both for neutrinos and antineu- 
trinos. In the same figure the axial-vector interference 
term is displayed as well as the isovector component il- 
lustrating the importance of the latter in the antineutrino 
case. A remark on the collective nature is in order at this 
stage. One notices the suppression produced by the RPA 
which follows from the repulsive character of the residual 
interaction. It is more pronounced at small energies. The 
interference term which spreads over all the energy range 
is therefore less affected by RPA. The antineutrino cross 
section which peaks at low energies is instead very sen- 
sitive to the suppression by the collective RPA effects. 
Once integrated over the energy transfer the reduction 
at a typical energy £'„=700 MeV is somewhat larger for 
antineutrinos than for neutrinos. A comparison between 
the two is affected by this difference. This could lead to 
a source of uncertainty in the comparison as the RPA 
effects which depend on the residual interaction have not 
been tested in this momentum regime. However we will 
show below that it does not prevent the comparison to 
be significant for our purpose. 

We now turn to the generalized "quasielastic" cross 
section which includes the multinucleon contribution. 
The neutrino and antineutrino genuine and generalized 
"quasielastic" cross section are plotted in Figj5]both in 
the bare and RPA case. As was already discussed in 
Ref.[ll[ the agreement with the MiniBooNE experimen- 
tal neutrino data 0, [Toj is better when the np — nh com- 



ponent is added to the genuine QE cross section, whether 
in the free or in the RPA case. Our prediction for the 
generalized neutrino "quasielastic" cross section shows 
only a moderate sensitivity to the collective aspects. For 
antineutrinos instead the sensitivity to RPA is somewhat 
larger but it does not hide the important point that the 
relative importance of the 2p — 2h term is smaller for 
antineutrino. This is illustrated in Fig|5] which shows 
the ratio of the multinucleon component to the single 
nucleon one with and without RPA. In both cases, we 
find that the ratio for antineutrino is reduced as com- 
pared to the neutrino one by a factor 1.7 at E v =700 
MeV. In order to eliminate the uncertainties related to 
the neutrino energy reconstruction, we give in Table 1 
values of quasielastic and multinucleon cross sections on 
12 C, as well as their sum, averaged over the respective 
neutrino [l(| and antineutrino [l7| MiniBooNE fluxes, so 
as to provide quantities more in touch with an experi- 
mental analysis. We gives these quantities both in RPA 
and in the free case and various situations are possible. 
For instance, if our RPA description holds, the ratio of 
the generalized "quasielastic" cross section, which is the 
measured cross section, to a theoretical free Fermi gas 
model is 1.22 for neutrinos and 0.99 for antineutrinos, 
significantly lower. If the extreme case where RPA ef- 
fects are totally absent the corresponding ratios are 1.37 
for neutrinos and 1.25 for antineutrinos. In all cases the 
antineutrino numbers are smaller and the difference may 
be detectable, which offers an experimental test. For 
neutrinos the fit to the quasielastic data in a relativistic 
Fermi gas description req uired an appreciable increase of 
the axial cut-off mass [10( ■ For antineutrinos the increase 
needed to account for the data in the same model should 
be smaller since the relative role of multinucleon ejection 
is reduced. This difference offers a possible way to shed 
light on the origin of the anomaly. Of course the differ- 
ence which occurs owing to the fact that the target is 
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QE np-nh QE+np-nh 


V 

QE np-nh QE+np-nh 


bare 


7.46 2.77 10.23 


2.09 0.52 2.61 


RPA 


6.40 2.73 9.13 


1.60 0.47 2.07 



TABLE I: MiniBooNE flux-integrated CC u^- l2 C and ^ M - 12 C 
total cross sections per neutron and per proton respectively 
in unit of 10~ 39 cm 2 . The experimental CCQE f M - 12 C value 
measured by MiniBooNE is 9.429 x 10" 39 cm 2 with a total 
normalization error of 10.7 % PH. 



not pure carbon but CH 2 affects exclusively antineutrino 
cross section reducing trivially the importance of nuclear 
effects. It has to be taken into account to reach a sig- 
nificant comparison between neutrino and antineutrino 
cross sections. 

Finally we would like to comment on the absence of 
final state interactions in our evaluation. In addition 
to the absorption of a produced pion it also ignores the 
possibility for an ejected nucleon to interact with the nu- 
cleus emitting another nucleon which leads to a final state 
with two nucleon ejected, the same type of final state as 
was discussed in this work. For the final sate interac- 
tion effect all responses are concerned and not only the 
spin isospin ones. This has been taken into account by 
Benhar and Meloni [l8| using a spectral function which 
describes the single particle dynamics in an interacting 
system, deduced from experimental (e, e'p) data. Their 
conclusion is that this inclusion cannot explain the en- 
hancement of the measured neutrino quasielastic cross 
section. Although the final sates are the same, (2p — 2h), 
we are dealing in the present work to a different type of 
correlations. These are ground state correlations, mostly 
tensor ones which affect only the spin isospin responses, 



producing an enhancement of the spin isospin sum rule 
which is reflected in an increase of the "quasielastic" cross 
section. The rest is due to the A excitation or the inter- 
ference between the two. These also obviously belong 
exclusively to the spin isospin sector, which is the ba- 
sis of our test through a comparison between neutrinos 
and antineutrinos. A fully realistic calculation should in- 
clude both the spectral function effect and the ground 
state correlations or A ones. 

In summary we have studied the quasielastic neu- 
trino and antineutrino cross section in the case of the 
MiniBooNE experiment where multinucleon ejection is 
not distinguishable from single nucleon production. The 
"quasielastic" cross section thus defined contains a cer- 
tain proportion of 2p — 2h and 3p — ih excitations. This 
proportion is large for neutrinos, which may be the inter- 
pretation of the increase in the axial cut-off mass needed 
to describe the data in the relativistic Fermi gas. For an- 
tineutrinos we predict a smaller role of the 2p — 2h com- 
ponent. The reason is that the vector- axial interference 
term produces a suppression of the spin isospin response 
contribution to the cross section leaving a larger role for 
the isovector response which is not affected be the 2p — 2h 
component. An experimental confirmation of this differ- 
ence would signal the fact that the excess cross section 
belongs to the spin-isospin channel, thus displaying the 
role played by the tensor correlations in neutrino nucleus 
interactions. 

The antineutrino mode which is actively investigated 
for the general problems of neutrino oscillations and CP 
violation is also of great relevance for the understanding 
of the neutrino nucleus interactions. 

We thank D. Davesne, T. Ericson, G.T. Garvey and 
T. Katori for useful discussions. 
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